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EFFECTS OF CAVITATION ON UNDERWATER SHOCK LOADING - 
PLANE PROBLEM, Part 2 

1. Introduction 

Development of a finite element fluid code for solving 
plane problems involving shock-induced cavitation was described 
in Ref. 1. The program uses displacement potential as the 
dependent variable as described in Ref. 2. A structural program 
(Ref. 3) providing a plane model for a ring-stiffened cylin- 
drical shell has recently been completed. In the following 
sections, coupling of the structural and fluid programs is 
described and results thus far are summarized. 



2. Structural Program 

2.1 Program Features . The structural program (STRUK5). was 
developed by Jack T. Waller, a thesis student of mine, and 
is described in Ref. 3. A shell structure having specified 
radius, area per unit axial length, radius of gyration, mater- 
ial density, and modulus of elasticity is modeled. Trigono- 
metric series are used to represent the radial and tangential 
components of shell displacement. Taking account of symmetry, 
the half-circumference of the shell is divided into N equal 
segments. The fluid mesh is matched with the structural 
model, providing an interface node at each point of division. 
Using fluid pressures calculated at these nodes by the fluid 
program, the coefficients for trigonometric series represen- 
tation of the interface pressure are found by STRUK5 . 
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For each trigonometric order m there is a term c^cosm0 
for interface pressure, a term a^^cosmS for radial displace- 
ment and a term b^^^sinme for tangential displacement. Because 
the structure is rotationally symmetric, there is no coupling 
between coefficients with different orders m. For a given 
m, differential equations governing a^ and b^^ are coupled. 

Time integration is performed using an explicit central 
difference algorithm. A chosen time step, which may be 
suitable with a particular fluid mesh, may result in instabi- 
lity for structural modes with large m. Accordingly the 
program limits the number of structural modes retained. For 
this purpose the extensional modes are considered separately 
from the flexural modes . 

2.2 Nonstructural Mass . The shell structure of a submarine 
typically has a mass no greater than 20 percent of the dis- 
placed water. The remaining mass is connected to the shell 
through relatively flexible structure and/or isolation mounting 
Although the early-time response to shock loading is little 
affected by this additional internal mass, its omission from 
the structural model will be significant at later times. 

Since cavitation effects on shock response may occur at late 
times, it was deemed important to include the additional mass. 

A rigid body displacement of the shell can be described 
by trigonometric terms of order m=l. For modeling purposes 
the spring forces between shell and contents are distributed 
in corresponding fashion. Thus there is a set of radial 
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forces varying as cos0 and a set of tangential forces varying 
as sinQ . This produces equations of motion in which the 
displacement of the nonstructural mass couples only with 
the shell displacements described by a^ and b^. 

Modification of the equations of motion to include this 
nonstructural mass ("contents" mass) was accomplished in 
the following way. The notation of Waller (Ref. 3) is used 
with the following additions. 

u^ = contents displacement (in 8 = 180° direction) 

m = contents mass 
c 

= common stiffness of radial and tangential spring sets 
The increments to elastic potential energy and kinetic energ-y, 
AU and AT, are given by 



AU = ?sk^[(a^-u^)^ + (b^-u^)^] 



AT = ijm u ^ 
c c 



( 1 ) 

( 2 ) 



When these additions are included the motion equations for 
a^ and b^ each contain both a^ and b^. Algebraic manipulation 
allows the equations to be written in a form suitable for 

explicit time integration as 

2 , c, 

(l+2z)a,+ ^( 1+z ) (a^ +b, )+— ^^[ ( 1+z ) a., -zb, - ( l + 2z )u. ]= (1+z )— ^ (3) 
1 11 iTpAR 11 c pA 

2 z c 

(l + 2z)b,+^(l+z) (a,+b, )+-^[-za, + (l+i)b, + (l+2z)u^] = --^ (4) 

^ X -L TT pAi\ XI C 

The additional equation of motion for the contents mass is 
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( 5 ) 




The contents mass m^ is determined within the program from 
the stipulation that shell plus contents be neutrally buoyant. 
Instead of specifying directly, the fixed-base frequency 
f of the contents is specified and k is found from 



2.3 Curvature Correction . When the just-described process 
for including nonstructural mass in the model was first 
implemented it was discovered that effects of shell rotatory 
inertia appeared in the equation for rigid body motion of 
the shell (obtained by putting b^=-a^ in the m=l equations). 
Since this behavior is clearly incorrect, the equations were 
carefully re-examined to discover the cause. Two related 
errors were discovered. Both appear first in my own notes 
and are clearly my responsibility, not Waller’s. 

The first mistake is that the expression used for 
rotational velocity u of the shell is incomplete. The 
correct form is 

^ " R 36 “ R 

The second term is omitted in Ref. 3. The other mistake is 
in the expression for shell bending curvature k . The correct 
form is 



c 



k = f ^/(SiT^m ) 
c c c 



( 6 ) 




( 8 ) 
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In Ref. 3 the second term in parentheses is +w. The latter 
form is correct only if the ring deformation is inextensional . 

i 

The reason that these errors were not discovered during 
initial testing of Waller's program is that they have negli- 
gible effects on results. Nevertheless, the corresponding 
changes have been made in the structural program. The new 
version of the program, which also includes the nonstructural 
mass as detailed in Section 2.2, is called STRUK6 . 

Revisions to the equations of Ref. 3 resulting from the 
corrections to rotational velocity (Eq. 7) and curvature 
(Eq. 8) are given in Appendix A. 

3 . Fluid Program 

Principal features of the fluid program were described 
in Ref. 1. Additional details are given in Ref. 4, including 
discussion of the staggered solution procedure which is used 
to couple fluid and structure for a time marching solution. 

4. Results 

Two kinds of results have been obtained from program 
tests conducted to date. The first kind give no new infor- 
mation concerning cavitation effects, but simply demonstrate 
that the coupled response of the structure is in agreement 
with known analytic results. These data, which afford 
supporting evidence concerning the validity of the modeling, 
are summarized in Appendix B. 
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Results of the second kind concern effects of cavitation 
upon the severity of shock induced structural loads. For 
the cases studied to date results are summarized in Ref. 4 
and quoted below. 



Examination of results from a large num- 
ber of cases studied up to this time reveals 
that the ratio of shock pressure increment to 
hydrostatic pressure (p^/p^) must exceed three 
times the ratio of decay length to structure 
radius (L/R) in order to induce cavitation. 
Cavitation begins about one decay length in 
front of the structure at about the time when 
the initial scattered wavefront reaches this 
location . 

Extreme values of bending moments and 
circumferential compressive force in the ring 
have been found to occur before the shock 
front has completely traversed the structure. 
Since cavitation has not yet begun during 
this interval, the extreme structural loads 
are not affected by cavitation. There is, to 
be sure , a shock pulse generated by cavity 
closure. The initial magnitude of this pulse 
is of the same order as hydrostatic pressure. 
After travelling to the structure the pulse 
is significantly attenuated and the resulting 
loads are not significant. 



5 . Conclusions 

Cases investigated thus far have shown no instance in 
which the occurrence of cavitation increases the maximum 
stresses in the structure. A wider range of encounter 
parameters needs to be investigated before it is possible 
to assert with confidence that cavitation induced by structure 
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response to underwater shock will not increase the likeli- 
hood of structural damage or collapse. 

Considering the related phenomenon of bulk cavitation, 
evidence to date does not rule out the possibility that it 
may intensify shock-induced structural loading. The finite 
element method should be used to investigate effects of 
bulk cavitation on structural loads. 



I 
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APPENDIX A 



Corrections for Equations of Reference 3 

Corrected forms for equations of Ref. 3 based on expression: 
for oj and < given in Eqs . 7 and 8 of this report appear below. 
Equation numbers are those of Ref. 3. 




(7) 

(9) 



dT = |p(ARd9)(w^ + v^) + ^IRde(^- 



5T = 7TpAR^(a^6a^ + &^5b^) + ^I(na^ + 



(1 + n^z)a^ t nzb; + (g(l + n^z)a^ + (f) n(l t n^z)b^ = ^ 

nza^ + (1 + z)b^ + (g'n(l + n^z)a^ + (f}'n^(l + z)b^ = 0 



\ 



(14) 



^ S„ * (f) “(n= - l)=za„ = '-i 

p2 n^pj '.R^ n^P^ pA 



b„ = - a„ /n 

"(F) '’(F) 



(17) 
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Equations 14 each contain both a and b and are thus unsuit- 
^ n n 

able for use with an explicit time integration algorithm. 

Some algebraic manipulation produces the following forms 
which are used in STRUK6 . 



(I + 2 + n^z)3„ t t n^z + (n^ - l)^z]a„ t (^'(1 + z)nb^ • (1 + z)^ 

(I t Z t „2z)b^ t (|“(1 . z + 2n"z)na„ + (||"(i t n“z)n\ • -r>2^ 
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APPENDIX B 



Modal Frequencies of a Submerged Cylinder 

It is well known that late time response motions of 
submerged structures can be accurately forecast by treating 
the surrounding fluid as incompressible (Ref. 5). This 
permits calculation of the frequencies for modes involving 
either rigid body translation or shell flexure. Comparison 
between such calculated results and observed frequencies 
of the shock response motion of the finite element model 
provides a partial check on the accuracy of the modeling and 
the implementing programs. 

Figure 1 shows the rigid body displacement a^ of the 
shell and the displacement u^ of the contents mass as a 
function of time in response to a shock with step pressure 
rise of 4MPa and decay time of 6ms. In view of the rapid 
shock pulse decay, the motion following the first half 
cycle may be considered free vibration. The succeeding 
cycle of motion has a period of 77ms or a frequency of 13.0 Hz. 
For the shell translation mode the added fluid mass is equal 
to the displaced mass. Taking this into account, the fixed 
base frequency f^ = 10.0 Hz for the contents mass gives a 
corrected value of 13.1 Hz in good agreement with the 
observed result. 

Figure 2 shows results for a separate computer solution 
with f^ = 5.0 Hz. In this instance the observed period is 
152ms, corresponding to a frequency of 6.6 Hz. This compares 
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Fig. 1. Shell and contents motion, = 10 Hz. 




Fig. 2. Shell and contents motion, = 5 Hz. 
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closely with a calculated value of 6.5 Hz. 

For the flexural modes corresponding to the trigonometric 
orders n = 2, 3, and 4, the corresponding added mass, based on 
maximum radial displacement amplitude, is (displaced mass)/n. 
Table 1 compares theoretic frequencies with finite element 
method (FEM) results. It is apparent that agreement is 
satisfactory . 



Table 1. Theoretic and Observed Modal Frequencies 



Order n 


Frequeni 

Theory 


::y (Hz) 

FEM 


2 


7.1 


7.3 


3 


23 


22 


4 


50 


48 
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